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Abstract 

The staggered magnetization of the Heisenberg antiferromagnet in 
two dimensions can be systematically approximated by a 1/N expan- 
sion. Cancellation between self energy diagrams leads to a Luttinger- 
like theorem for the ground state. We prove (for a smooth enough self 
energy) that the long range order of mean field theory (iV=oo) survives 
corrections to all orders of 1/N. Divergences of this series provides a 
new route to the disordered phases of quantum antiferromagnets. 
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In the study of quantum phase transitions, the order-to-disorder transi- 
tions of the Heisenberg antiferromagnet at zero temperature are particularly 
interesting. There is also hope that understanding such transitions may 
provide insight into the electronic correlations of lanthanum cuprates where 
under low doping, antiferromagnetism is replaced by superconductivity. 

The ground state of the Heisenberg antiferromagnet in two dimensions 
can have either long range order or be disordered by quantum fluctuations 
U]. Disordering can be induced by frustrating longer range interactions, 
or perhaps by slight hole doping as suggested by the phase diagram of lan- 
thanum cuprates. In either case, there are theoretical difficulties in describing 
the transition itself. While semiclassical spin wave theory works well deep 
in the ordered phase, it assumes spontaneously broken symmetry, and fails 
when the staggered magnetization becomes small. The continuum approx- 
imation of the (2+1) dimensional non linear sigma model relies on perfect 
short range antiferromagnetic correlations. Near the transition however, the 
short range correlations deteriorate considerably. This complicate matters, 
since one needs to include field discontinuities (e.g. hedgehogs), and consider 
interference effects between their Berry phases ||. 

The Schwinger Bosons (SB) large-N expansion || is a rotationally sym- 
metric formulation, which in principle can treat both sides of the transition 
[|]]. The mean field theory (MFT) describes the excitations as a free Bose 
gas of iV decoupled flavors. Bose condensation in this system is equivalent 
to long range spin order ||. However, MFT is strictly valid only at N = oo, 
while the physically interesting system is at N=2. A connection between the 
two limits requires an understanding of the 1/N series. Higher order cor- 
rections involve interactions between SB which enforce the local constraints. 
However, finite N corrections to the staggered magnetization have not yet 
been evaluated. Before this could be done, it is necessary to place the 1/N 
expansion on firmer footing, i.e. to show that the higher order terms yield 
finite and sensible results, which do not immediately destroy the mean field 
ground state. Currently, we do not know whether the long range order found 
in the MFT survives for any N < oo. 

This paper specifically addresses this concern. We prove a theorem which 
establishes the 1/N expansion as a consistent approach for the ground state 
of finite N systems, starting from the MFT. Under a condition that the self 
energy is sufficiently smooth at the ordering wavevectors, we prove that if 
there is long range order in the MFT, the spontaneous staggered magnetiza- 
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tion does not vanish to all orders of the 1/N expansion. The proof uses a 
cancellation between self energy diagrams and their tadpole counterparts, a 
feature special to the 1/N expansion. 

The result is reminiscent of (but not equivalent to) Hugenholtz and Pines' 
self energy condition for Bose condensed liquids 0. It is closer in spirit 
to Luttinger's theorem for Fermi liquids 0. The spontaneous staggered 
magnetization is analogous to Luttinger's Fermi surface. (Both appear as a 
discontinuity in the occupation number.) In Luttinger's theorem, under a 
similar condition on the self energy, the Fermi surface discontinuity survives 
at each order in perturbation theory. Pushing this analogy further, we shall 
propose that the vanishing of the staggered magnetization at finite N may 
formally resemble one of the known Fermi surface instabilities. 

For simplicity we discuss the nearest neighbor SU(N) Heisenberg anti- 
ferromagnet. The proof actually utilizes only general features of this model, 
and thus it is readily extendable to more general hamiltonians. The spins are 
represented by N SB per site a\_ mmax . . . a\ mmax where m max = (N - l)/2, 
and the Hamiltonian is given by 

H = — — ^{a\ m a) m )(a im ,a jrril ) - hJ2 ma\ m a im , (1) 

(i,j) are nearest neighbor bonds on the square latice. Summation over 
repeated indices is implicit, unless specified otherwise. The Hilbert space 
is constrained by the fixed SB number a\ m ai m = Ns, at each site, h 
is an infinitesimal ordering field. The generators of SU(N) are given by 
a>lm a im', where we use conjugated representations on opposite sublattices. 
For N=2, Eq. (|1|) is equivalent (up to a constant) to the Heisenberg model 

n = jj2 {i , j) s l -s 3 -hE t (-i) i S! !. 

Following the standard procedure p3] , the partition function can be writ- 
ten as a coherent states path integral. Hence one introduces real local fields 
Aj to impose the constraints, and Hubbard Stratonovich fields Qij to decouple 
the quartic interactions. 

After integrating out the SB field, we are left with 

Z{h) = fv(XQ) exp[-NS(\,Q,h)] (2) 

The explicit expression for the action S can be found in Ref. ||. Following 
|J, we consider the case of zero temperature and large, but finite volume. 
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The staggered magnetization is given by 

M = lim lim Tr mG(k, m) (3) 

where G is the full Green function of SB. It is given in the Nambu notation 
as a 2 x 2 matrix with normal and anomalous components G nn r(k,m) = 
(A n km^ n >k m ) > where A km = (a km , a_ km ). (3 is the inverse temperature and V 
is the volume (number of lattice sites). We denote k = (k, u), where k and 
uj are lattice momentum and Matsubara frequency respectively. Tr includes 
a trace over k, m and the Nambu indices. G is evaluated by summing all 
one-particle diagrams generated by the large-N expansion of Eq. (|2]). 

Let us briefly review the mean field results which were derived previously 
[||, |[. At large N, (0) is dominated by the saddle point A, Q, and G is 
approximated by the mean field Green function Go, 

G ^ m )={ 4 Q 7k A + ^-Z^J (4) 
where 7k = ^(cosk x + cosk y ). The poles of Go are at the SB frequencies 



\| A ' + ^ h ( mm ™ ~ m ) + 2 ^ - 7 k) (5) 



where c = V8Q and = c 1 y\ 2 — 16Q 2 — 2h\m max . uj^ is minimized at 
the two points k c = (0, if). At those momenta, for m = m max , the excitation 
gap is cA h . Solving the mean field equations yields A h = y/2/(NV(s — 
0.1966 . . .)). In the thermodynamic limit, the SB with m = m max and k = 
k c contribute macroscopically to the momentum sum, i.e. they undergo 
Bose condensation. This condensate is the only term, which survives the 
cancellation between positive and negative m's in (Q), yielding the mean 
field staggered magnetization 

M„^^^( S -0, 96 ,..) (6 , 

which for N = 2 agrees with spin wave theory 0. 

The higher order 1/N corrections to G are described by diagrams which 
include lines for Gq interacting via propagators D (defined later) which are 
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depicted as wiggly lines. A diagram which involves L loops (traces of prod- 
ucts of Go) an d P propagators, is of order {1/N) P where p = P — L. One 
must exclude all diagrams which include the segments shown in Fig. 1. As 
shown in 0, §|, this leads to the fulfillment of the SB constraints at each 
order oil/N separately. 

As in MFT, a nonzero staggered magnetization is related to the divergence 
of the number of SB with m = m max at k = k c . On the other hand, strictly 
at h = 0, MFT is SU(N) rotationally invariant, so the Bose condensation 
is equally shared among the different m flavors and the gap becomes Ao = 
NAhjLo- Henceforth we shall set h = 0, and have exact degeneracy between 
the different flavors (jj^ lfn . Thus, long range order is associated with Bose 
condensation of all flavors at k c . 

The self energy £(&) (also a 2 x 2 matrix) is related to G by the Dyson 
equation G -1 = Gq 1 — E. In order to proceed we must make an impor- 
tant assumption on the smoothness of the self energy near the condensate 
momenta: 

lim \Z(k) - Z(k c )\ = O (\k - k c \ 2 - s ) , 6<1 (7) 

k^k c \ / 

where k c = (k c ,0). S should exhibit rotational symmetry about k c as a 
consequence of the asymptotic "Lorentz invariance" of Go near k c . (The SB 
dispersion vanishes linearly at k c .) We have verified that the leading order 
self energy is smooth at k c (i.e. obeys (0) with 5 = 0). We argue that the 
smoothness assumption is plausible for models which have no pathology in 
the density of low excitations. For such models, the integrations in S are 
uniformly convergent for all external momenta. 

However, we have not proven Eq. (|7|) to all orders in 1/N, and we must 
regard it as an assumption; one which requires a separate justification for 
any particular model. 

The number of SB with momentum k c is = (2j3)~ l J2u TrG(k c , u), 
where Tr traces over Nambu indices. This number diverges as n^ c ~ V if 
det[G _1 (A; c )] ~ V~ 2 . We use the Dyson equation and MFT relation A ~ 
V^ 1 to state that 

A' = E u (k c ) - 7kc S 12 (A; c ) = 0(^ 2 ) ~> M ± (8) 

i.e. if the quantity A' vanishes rapidly enough in the thermodynamic limit, 
the ground state has long range order. It may be shown, that A'(0) = A' (7?). 
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Theorem: Under condition (^), Eq. (^j) holds to all orders in the 1/N 
series. 

Proof: The self energy is decomposed into two parts 

E = E + S* ad (9) 

where S* ad is the single tadpole diagrams (see Fig. 2), and £ are all the 
remaining diagrams. Although £ and E tad are expected to be of 0(1) sepa- 
rately, we shall show that at k = k c the 0(1) contributions precisely cancel 
in (ID leaving us with terms of 0{V~ 2 ). Note that in contrast to perturba- 
tion theory, the first and second terms in Fig. 2. have a different number 
of vertices, but are of the same order in 1/N. This enables the cancellation 
mechanism function at each order separately. 

The rest of this discussion contains unavoidable technical details. The 
set of auxiliary fields is denoted by (Xj, 3tQ s ,j, ^Qd,ji ^Qs,j, ^Qd,j)- \ cou- 
ples to the local Boson density and dt(^s)Q s ^d),j couple to the bilinear forms 
Y J e=e x ,e y r n S e {d) [a)a) +e + Ha^+J, where rj^ = r) s By = -qj x = 1 and rf ey = -1. 
We define 2x2 vertices v a which connect between a field a and two Go's. 
Thus a zero momentum field a is coupled to the form Y^k^v^Ak, where 

= ^k' 3 = cr x (cosk x ± cosky), and -O^ 5 = ia y (cosk x ± cosk y ). Using 

v a , we can explicitly write £* ad (/c c ) (see Fig. 2) as 

Y. tad (k c ) = 2NvlD a ' a '{0)J2 r ^[^G (k)R(k)G (k) 

k 

R = £ + £G£ (10) 

where Y,k = y Sk / f^r- It ma y be seen, that only a, a' = 1, 2 give nonvan- 
ishing contribution to this formula. 

£(&) may be expanded as £(&) = T, a {k)u^, where summation over a runs 
from 1 to 3, u^ 2 = v^ 2 and u k = o z . The coefficients of expansion satisfy 
the relation E Q (0, 0) = £ Q (7r, 0) and T,^(kc) = 0. The same expansion with 
u a holds for R(k). A' can now be written as A' = f a H a (k c ), where f 1 = i/2, 
f 2 = -2 and f a = 0, a > 2. 

The propogator in the 1/N expansion is given by the matrix D = 
j, (n - n) _1 , where 

4 

Uo' a (q) = 5 a>a >(l - <W)j 

U^'(g) = 2^Trk k+q Go(fc + g)^; qik G (A ; )] (11) 



6 



Using the relation NDU = — 1 + NDUq, we find that at k = k c , S(fc c ) is 
cancelled on the right hand side of Eq. (Q) and we obtain 



Y^ a (k c 



AN 

- [Z{k c )G(k c )X(k c )] a + —D a > 2 (0)R 2 (k c 



x 



1 - G (k)E(k)J a „ \l - G (fc c )£(fc c ), 



(12) 



where we have used the fact that is independent of momentum. 

The Bose condenstation of Go(k c ) gives rise to the divergence of IT(g = 0). 
Extracting the volume divergences in (|n|) yields Tl a > a '(0) = f a f a '(aV 2 + 
bV) + P a,a where a, b, and P a,a are independent on volume. We see that 
both order V 2 and order V factorize, reflecting the emergence of a discon- 
nected part in the correlation function due to Bose condensation. Denoting 
P = n — P and inverting the polarization matrix we obtain the propagator 
to order V' 2 : 



ND a > a '(0) = (P 



-l\a,a' 



/7(p- 



1^7,7' fl' 



aV 2 [/7(p-i)7,7' f-y'Y 

(13) 

We now expand S in a power series of 1/N: S = J2 p N P S^. We shall 
prove Eq. (§) by induction. Assume that Eq. (|8|) holds for SW,p < p. 
We take on the left hand side of Eq. ([12]), and multiply both sides 

by the vector f a . Using the Dyson equation for G, one can show, that 
/ Q [SGS] Q is proportional to A', which, however, should be calculated with 
YP-P. Therefore, this term yields 0(V~ 2 ). Then, using ([13]), the terms of 
0(1) in D get cancelled by multiplying them on the left by f a , leaving us 
with an overall factor of 0(V~ 2 ). We must still show that the second factor 
in the third term of ( ]T2| ) is not divergent. Since the summand diverges as 
{k — k c )~ 2 , (2 powers of the phase space minus four powers from Go), the 
momentum sum will converge if the self energy obeys condition (0). Thus, 
we have shown that Eq. (§) holds to all orders in 1/N. Q.E.D. 

We note that it is crucial for the cancellation, described above, that the 
constraint has a local character and enforced by a fluctuating field. Indeed, 
this cancellation does not take place, if the constraint is imposed only on 
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average by a static chemical potential. On the other hand, our proof can 
be readily extended to different spin models with constrained Hilbert spaces. 
In particular it applies to the t'-J model, a semiclassical approximation to 
holes in the quantum antiferromagnet ||. Also, a simpler version of this 
theorem applies to the long range order in Resonating Valence Bonds states 
[PHI , using a large-N expansion of the Gutzwiller projection ||. 

In practical terms, this theorem sets the foundation for investigating the 
disordering transition using the 1/N expansion of the self energy. We can 
propose two scenarios for the disordering mechanism at finite N: (i) Coupling 
of spins to soft charge fluctuations (holes) can give rise to violation of (|7[), i.e. 
a breakdown of our theorem and a destruction of long range order. This sce- 
nario is analogous to the one dimensional Luttinger model where the Fermi 
surface discontinuity vanishes due to the large density of low excitations, (ii) 
The divergence of V 2 Y, p N~ P T,^ — > oo may be detected in a partial resum- 
mation scheme. (Tadpole countertems must be properly included, as shown 
above.) A divergence for example in nested diagrams, formally resembles the 
Cooper channel (superconductivity) instability in a Fermi liquid. 

In summary, we have analyzed the corrections to the mean field ground 
state staggered magnetization of the two dimensional antiferromagnetic Hei- 
senberg model. We found an important cancellation mechanism between self 
energy diagrams. This establishes that the 1/N expansion for the order pa- 
rameter is a consistent asymptotic approach for finite iV models. It is similar 
to perturbation theory about a non interacting Fermi surface. We argue 
that the quantum disordering transition may be detected as a breakdown of 
the assumptions of this theorem, or a divergence in the 1/N series. These 
possibilities are worth further investigations. 
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Figure Captions 

Fig. 1: Forbidden segments in 1/N expansion diagrams. Solid lines repre- 
sent mean field Green functions Go and wavy lines - propagators of 
auxiliary fields D (see text). 

Fig. 2: Diagrammatic representation of Eqs. (P, 0) for the self energy. E 
are all diagrams except the single tadpole diagrams. The cancellation 
between S and the tadpole diagrams allows the staggered magnetiza- 
tion to survive finite N corrections. 
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